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Molecular line intensities as measures of cloud masses - 
I. Sensitivity of CO emissions to physical parameter 
variations 
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ABSTRACT 

A reliable estimate of the molecular gas content in galaxies plays a crucial role in 
determining their dynamical and star- forming properties. However, H2, the dominant 
molecular species, is difficult to observe directly, particularly in the regions where 
most molecular gas is thought to reside. Its mass is therefore commonly inferred by 
assuming a direct proportionality with the integrated intensity of the 12 CO (J = 1 — ► 0) 
emission line, using a CO-to-Eb conversion factor, X . Although a canonical value for X 
is used extensively in such estimates, there is increasing evidence, both theoretical and 
observational, that the conversion factor may vary by over an order of magnitude under 
conditions different to those of the local neighbourhood. In an effort to understand 
the influence of changing environmental conditions on the conversion factor, we derive 
theoretical estimates of X for a wide range of physical parameters using a photon- 
dominated region (PDR) time-dependent chemical model, benchmarking key results 
against those of an independent PDR code to ensure reliability. Based on these results, 
the sensitivity of the X factor to change in each physical parameter is interpreted 
in terms of the chemistry and physical processes within the cloud. In addition to 
confirming previous observationally derived trends, we find that the time-dependence 
of the chemistry, often neglected in such models, has a considerable influence on the 
value of the conversion factor. 
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1 INTRODUCTION 

Accurately measuring the mass of molecular gas in galaxies 
is of great importance in determining the dynamical proper- 
ties and star formation efficiency that are critical to models 
of galaxy evolution. Such measurements also play a vital role 
in understanding the properties of the interstellar medium 
in extragalactic environments. Molecular hydrogen (H2) is 
the dominant molecular species in the Universe and governs 
the dynamical and chemical evolution of molecular clouds. 
However, its lack of a dipole moment means that it can only 
change ro-vibrational state through weak quadrupolar tran- 
sitions that are only rarely detectable. Alternative means of 
deriving the mass of molecular clouds have therefore been 
developed over the years. Many of these techniques make 
use of a tracer species (a molecule or dust) that is easier 
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to detect and infer the molecular hydrogen content by some 
proportionality relation. 

One of the most commonly used tracers of molecular gas 
is carbon monoxide (CO). After H2, it is the most abundant 
molecule in molecular clouds by several orders of magnitude 
and has been observed in all regions where molecular gas 
is thought to be present. However, unlike molecular hydro- 
gen, CO readily emits under the conditions typically found 
in molecular clouds and is easily observed at radio and sub- 
millimetre wavelengths with ground-based facilities. 

The column density of molecular hydrogen along a line 
of sight is often related to the integrated intensity of the 
12 CO(J = 1 — > 0) rotational transition line using the CO- 
to-H2 conversion factor: 

A(H 2 ) 



X = 



J T A (CO)di 



[cm 2 (Kkms 1 )~ 



(1) 



where A(H2) is the column density of H2 and Ta(CO) is the 
antenna temperature of the CO(1-0) line (see, e.g., reviews 
by Maloney 1990; Combes 1991; Young & Scoville 1991). 
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Observational estimates of X have been derived for 
nearby regions where other methods for determining molec- 
ular cloud masses are available. Values of X were first de- 
termined using extinction measurements obtained from star 
counts to derive N(H.2), assuming a constant gas-to-dust 
relationship, and comparing with observed CO luminosities 
(e.g. Dickman 1978; Frerking, Langer & Wilson 1982). Scov- 
ille et al. (1987) and Solomon et al. (1987) derived values 
of X for a large sample of resolved giant molecular clouds 
with assigned kinematic distances by calculating their virial 
masses from CO linewidths. The same approach has re- 
cently been applied to clouds within nearby galaxies (e.g. 
Rosolowsky et al. 2003; Rosolowsky & Blitz 2005). It should 
be noted that the values derived using this method are based 
on the assumption that the clouds are virialized. Maps of the 
diffuse 7-ray emission produced by interactions between cos- 
mic rays and hydrogen nuclei have been compared to H I and 
CO survey data for the Galactic plane to calculate the con- 
tribution of H2 nuclei to the 7-ray emission and so calculate 
X values along many lines of sight (Blitz et al. 1985; Bhat, 
Mayer & Wolfendale 1986; Richardson & Wolfendale 1988; 
Strong et al. 1988; Bloemen 1989; Strong & Mattox 1996). 
The remarkable outcome of these studies is that the value 
of X seems to be roughly constant for the Galactic molecu- 
lar clouds that have been observed, and a canonical value is 
often adopted: X w 2 x 10 20 cm" 2 (Kkms^ 1 )" 1 (Strong & 
Mattox 1996; Dame, Hartmann & Thaddeus 2001). There 
have been several attempts to explain this apparent invari- 
ance, citing ensembles of virialized clouds (Dickman, Snell 
& Schloerb 1986; Solomon et al. 1987) and dumpiness of 
the interstellar medium on scales of Ay < 2 mag (Taylor, 
Hartquist & Williams 1993). 

Despite the widespread use of this canonical value, there 
is strong evidence to suggest that X can vary significantly 
under conditions differing from the local environment. Even 
within the Milky Way, observationally determined values for 
the X factor are found to vary with Galactocentric radius 
(Digel, Bally & Thaddeus 1990; Sodroski et al. 1995; Dah- 
men et al. 1998) and, further afield, estimates for X in the 
Small and Large Magellanic Clouds are over an order of mag- 
nitude higher than the canonical value (Cohen et al. 1988; 
Rubio et al. 1991; Israel 1997). Conversion factors up to an 
order of magnitude lower than the canonical value have also 
been determined for a number of galaxy nuclei (Smith et 
al. 1991; Israel & Baas 2001; Israel, Tilanus & Baas 2006). 
Theoretical studies of the variation of the CO-to-Hb ratio 
have also suggested that X can be up to an order of magni- 
tude higher or lower than the canonical value under certain 
conditions (e.g. Sakamoto 1996, 1999). Taylor et al. (1993) 
used a photon-dominated region (PDR) chemical model to 
perform ab initio calculations of X and showed its varia- 
tion for small changes in typical cloud parameters. Work 
along similar lines has also been conducted by Le Bourlot, 
Pineau des Forets & Roueff (1993a), Wolfire, Hollenbach & 
Tielens (1993) and Kaufman et al. (1999). These previous 
studies provide the foundation and motivation for this cur- 
rent project. 

Since the majority of CO emission occurs in the outer 
regions of molecular clouds, where the interstellar radiation 
field falls on the cloud surface, PDR models are most suit- 
able to the theoretical study of the X factor. PDRs arise 
wherever FUV radiation (6 < hv < 13.6 eV) is incident upon 



Table 1. Elemental abundances used in the UCL_PDR code (rela- 
tive to total hydrogen nuclei). 
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7.50x10" 


-2 


C 


1.42x10- 
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3.19x10- 
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Si 
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1.43x10" 


-6 


CI 


1.10x10" 




Ca 


5.72x10- 


-10 


Fc 


3.60x10- 











the outer edge of a molecular cloud. The radiation drives 
photochemistry within the cloud and governs the chemical 
and thermal structure near its surface. 

In this paper we examine the influence of physical and 
chemical conditions on X and attempt to explain its depen- 
dence on important physical parameters, such as the density, 
age, and metallicity of the cloud. In Section 2 we describe the 
PDR models used to calculate and independently verify the- 
oretical estimates of the CO-to-H2 ratio for given physical 
conditions. Section 3 reviews the independent benchmark- 
ing results for important combinations of parameters when 
using the two PDR models and evaluates the reliability of 
subsequent, more detailed calculations. The effect of varying 
physical parameters is then discussed in depth in Section 4, 
before we summarize our conclusions in Section 5. 



2 PDR MODELS 

2.1 The UCL_PDR code 

We use the ucl_pdr time-dependent PDR code (Pa- 
padopoulos, Thi & Viti 2002; Bell et al. 2005) to calcu- 
late the chemistry, thermal balance and CO (1-0) emission 
strength within the outer region (Av < 10 mag) of a cloud 
to derive theoretical values of the CO-to-H2 ratio, X, for 
given physical parameters. 

The PDR is modelled as a one-dimensional semi-infinite 
slab, illuminated from one side. The chemistry and ther- 
mal balance are calculated self-consistently at each depth 
point into the slab and at each time-step, producing chem- 
ical abundances, emission line strengths and gas tempera- 
tures as a function of depth and time. 

We adopt a chemical network containing 128 species 
and over 1700 reactions, including ion-molecule, photoion- 
ization and photodissociation reactions. Freeze-out of atoms 
and molecules onto grains is neglected. The gas is assumed 
to be initially in atomic form, with the elemental abun- 
dances listed in Table 1 (all fractional abundances in this 
paper are quoted relative to total hydrogen nuclei). Elemen- 
tal abundances of all metals are assumed to scale linearly 
with metallicity (Z/Zq). The reaction rates are taken from 
the UMIST99 database (Le Teuff, Millar & Markwick 2000), 
with some modifications introduced as part of a recent PDR 
benchmarking effort (Rollig et al. 2006). We adopt the grain 
surface H2 formation rate of de Jong (1977): 

R = 3 x 10~ 18 \/Texp(-T/T000) [cm 3 s _1 ] (2) 

and assume that the rate scales linearly with metallicity. 

Extinction within the cloud is calculated assuming an 
average grain size of 0.1 [im, albedo of 0.7 and mean photon 
scattering by grains of g — 0.9 (g = (cos 9}). The dust-to-gas 
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Table 2. Standard parameter values adopted for each cloud 
model. 



Parameter Standard Value 





10 3 cm 3 


X 


1 Draine 


C 


1.3 x icr 17 s- 1 


z 




VtuTb 


1 kms 1 


t 


1 Gyr 



mass ratio is assumed to scale linearly with metallicity and 
takes a standard value of 10 -2 at Solar metallicity (Z@). 
The visual extinction (Av) and dust optical depth (t v ) are 
related by: 

Ay = 1.086 t v (3) 

The incident FUV radiation is characterised by the 
standard Draine field (Draine 1978), scaled by a free pa- 
rameter x- Attenuation of the radiation field strength with 
cloud depth is given by: 

X = Xoexp(-fc^v) (4) 

where k — 1.38 and xo is the unattenuated field strength at 
the surface of the cloud (in units of the Draine field). 

Self-shielding of Eh against photodissociation is cal- 
culated using the single line approximation of Federman, 
Glassgold & Kwan (1979). The self-shielding factors of van 
Dishoeck & Black (1988) are adopted for CO. Neutral car- 
bon shielding includes contributions from H2, CO and Cl, 
using the treatment of Kamp & Bertoldi (2000). 

The gas temperature is determined at each depth point 
using an iterative process to balance the total heating and 
cooling rates. Photoejection of electrons from dust grains 
(Tielens & Hollenbach 1985) and PAHs (Bakes & Tie- 
lens 1994), and FUV pumping and photodissociation of H2 
molecules (Hollenbach & McKee 1979) are the dominant 
heating mechanisms near the cloud surface, whilst cosmic 
ray heating (Tielens & Hollenbach 1985) and exothermic re- 
actions become significant deeper into the cloud. Heating by 
Cl photoionization (Kamp & van Zadelhoff 2001), H2 for- 
mation (Kamp & van Zadelhoff 2001), gas-grain collisions 
(Burke & Hollenbach 1983) and turbulence (Black 1987) are 
also included in the model. 

The gas is mainly cooled through emission from col- 
lisionally excited atoms and molecules and by interactions 
with the cooler dust grains. Emission from the [Ol], [Cl] 
and [C 11] fine structure lines and CO rotational lines (up to 
J = 11 — * 10) are calculated at each depth point in the code 
using the escape probability method of de Jong, Boland & 
Dalgarno (1980) and non-LTE level populations determined 
using the updated collisional rate coefficients from the recent 
PDR benchmarking (Rollig et al. 2006) . The contribution to 
the total cooling rate by H2 emission is calculated using data 
from Martin, Schwarz & Mandy (1996). Ha and Ol 6300A 
fine structure line emission are also considered. 

Variation in metallicity is accounted for by assuming 
that the dust and PAH photoelectric heating rates are di- 
rectly proportional to Z, in addition to the scaling of the 
elemental abundances, H2 formation rate, dust-to-gas mass 
ratio and grain number density. 



Table 3. Benchmark model parameters. When not specified, the 
standard parameter values listed in Table 2 were used. 



Model 


Parameter 


A 


Standard Values 


B 


7iu = 500 cm -3 


C 


X = 10 3 Draine 


D 


C = 1.3 x 10~ 14 s- 1 


E 


Z = 10- 2 X Zq 


F 


"turb = 3 kms -1 



The integrated intensity of the CO (1-0) line along the 
line of sight towards the cloud is then: 

I— — / A(z) dz [erg s -1 cm -2 sr -1 ] (5) 
271 J 

where A(z) is the CO(1-0) line emissivity at depth z into 
the cloud (ergs -1 cm -3 ) and the factor of 1/27T accounts for 
the fact that photons only emerge from the outer edge of 
the cloud, rather than over the full 47t steradians. 

The theoretical antenna temperature, integrated over 
the line frequency, is given by: 

The value of X is therefore calculated at each depth 
point into the PDR using equation 1. 

2.2 The Meudon code 

The Meudon PDR code self-consistently computes the 
steady-state chemical and thermal structure of the cloud, 
in addition to the emission strengths of the relevant cooling 
lines. A new version of the Meudon code is now available (for 
details, see Le Petit et al. 2006), which includes a detailed 
treatment of the ultraviolet radiative transfer and the exten- 
sion to finite plane-parallel geometry. However, the present 
benchmark calculations have been performed using an older 
version of the code (Abgrall et al. 1992; Le Bourlot et al. 
1993b; Flower et al. 1994), in which self-shielding of H 2 and 
CO are computed using the FGK approximation (Federman 
et al. 1979) and do not include the overlap between H2 and 
CO predissociating lines. This has been demonstrated to be 
a reasonable approximation by Lee et al. (1996). 



3 BENCHMARKING 

In order to assess the reliability of the model predictions 
in this work, the X factor data produced by the UCL_pdr 
code have been benchmarked against the model results of 
the Meudon PDR code to check that the same trends are 
found in both. The Meudon code treats the microphysical 
processes within the cloud in more detail than the ucl_pdr 
code, but is more computationally intensive and less suited 
to examining large regions of parameter space. Therefore, to 
test the agreement between the two codes, we have selected 
six key models (listed in Table 3) that are representative of 
the parameter space being examined. 

The same species set and reaction network have been 
used in both codes as part of the benchmarking process. The 
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Figure 1. Benchmark model comparisons. X factor as a function of Ay, as calculated by the UCL_PDR (solid line) and Meudon (dotted 
line) codes, is plotted for models A (standard parameters; panel a), B (njj = 500 cm -3 ; panel b), C (x = 10 3 Draine; panel c), D 
(£ = 10 3 ; panel d), E (Z = 10~ 2 ; panel e) and F («turb = 3 kms" 1 ; panel f). 



grain properties, extinction law and metallicity-dependence 
are also common to both codes, so that the results can be 
directly compared. Since the Meudon code computes steady- 
state chemistry, a cloud age of 1 Gyr has been used in the 
UCL_pdr code to approach chemical equilibrium. 

Fig. 1 shows the results of the benchmarking tests. It 
must be stressed that exact agreement between the two 
codes is not expected, since their treatment of the physi- 
cal processes is very different in some cases. Overall, agree- 
ment is good between the two codes, with X profile minima 
generally differing by less than a factor of 2 across the vi- 
sual extinction range considered here (see Table 4) . The test 
models probe the influence of each parameter considered in 
this paper and both codes predict the same trend for the 
X factor variation in each case. However, the Meudon code 
produces X factor depth profiles with minima that are con- 
sistently lower and deeper into the cloud than those calcu- 
lated by the UCL_pdr code. This can be attributed to the 
different treatments of the self-shielding used by the two 
codes in the attenuation of the CO photodissociation and 
Cl photoionization rates (see Sections 2.1 and 2.2 for de- 
tails). The Meudon models also consider the increased pho- 
todissociation rate of H2 molecules in rotationally excited 
states, an effect neglected in the UCL_PDR code. The differ- 
ences caused by the individual self-shielding treatments and 
the additional photodissociation of excited H2 included in 
the Meudon code are most apparent in benchmark model 
C (x ~ 10 3 Draine; panel c in Fig. 1) where the FUV ra- 
diation is strong and the gas temperatures are high. These 
differences should be preserved, since they reflect the choice 
of treatment for this physical process and give an indication 
of the uncertainty inherent in calculations of theoretical X 
factors. We regard the agreement between the codes as sat- 



Table 4. Minimum X values obtained for the six benchmark 
models, as calculated by the UCL_PDR and Meudon codes. 



*min [10 20 cm- 2 (Kkms- 1 )" 1 ] 
Model UCL.PDR Meudon 



A 


0.9 


0.7 


B 


1.5 


1.0 


C 


1.4 


0.7 


D 


0.6 


0.8 


E 


300.0 


190.0 


F 


0.5 


0.4 



isfactory and the exploration of parameter space reported 
below is made using the ucl_pdr code. 



4 MODEL RESULTS AND SENSITIVITY OF X 
TO PARAMETER VARIATIONS 

We now examine the variation of the X factor under differ- 
ent physical conditions. The influence of changing the num- 
ber density of hydrogen nuclei (tih), incident radiation field 
strength (%), cosmic ray ionization rate (£), cloud age (t), 
metallicity (Z) and turbulent velocity («turb) are considered 
individually. With the exception of the parameter being ex- 
amined, the standard values listed in Table 2 are adopted 
for each model. Plots of X as a function of visual extinction 
(Av) into the cloud are used to illustrate these trends for 
each parameter. 

The characteristic depth profile of the X factor has 
three main features. Near the outer edge of the cloud, where 
the H/H2 transition occurs and H2 is present without CO, 
X is large. At higher visual extinction, while all hydrogen 
is now in molecular form, CO self-shielding becomes suf- 
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ficiently strong to prevent photodissociation, allowing its 
abundance and emission strength to rise sharply. This causes 
X to drop, reaching a minimum value in the region where 
CO emission is strongest. As the CO line becomes opti- 
cally thick with increasing depth, local emission gradually 
declines and X begins to rise slowly again. The minimum 
of the depth profile therefore represents the value of X that 
would be observed for that particular environment, since it 
corresponds to the peak emission of the CO line (see Fig. 2), 
with little contribution from deeper within the cloud. Under 
most conditions, X profile minima tend to occur in the re- 
gion Av < 2 mag where freeze-out of atoms and molecules 
should not take place, therefore our decision to exclude it 
from the chemistry seems reasonable. 

The lower panel of Fig. 2 demonstrates that once the 
CO line has become optically thick, its integrated intensity 
reaches a constant value and is no longer sensitive to the 
monotonically increasing H2 column density deeper into the 
cloud. Therefore, in the case of a slab of uniform density, 
the X factor only provides a lower limit for 7V(fi2). How- 
ever, in reality, there is good evidence that at least some 
interstellar clouds are clumpy (see, for example, Peng et al. 
1998; Morata, Chart & Estalella 2005). In this case, the CO 
emission arises from an ensemble of clumps, each effectively 
optically thin. The total H2 column density is then recovered 
using the X factor method if clouds are clumpy on scales of 
Av ~ 1 mag, as has been suggested previously (e.g. Taylor 
et al. 1993; Wolfire et al. 1993). 

For our choice of standard physical parameters adopted 
in Fig. 2, the derived value of X would be about 1, rather 
than the canonical value for the Milky Way of 2, in units of 



10" 



10 20 cm" 2 (Kkms" 



We note, however, that we could 



obtain X = 2 with a variety of plausible parameter combina- 
tions, e.g., by setting the cloud number density to 400 cm~ 3 , 
or the cloud age to ~1 Myr, or the turbulent velocity to 
0.3 kms -1 . This degeneracy in reproducing the observed 
value therefore shows that X cannot be used to infer the 
physical conditions and, conversely, that using the canoni- 
cal value in any situation without further consideration can 
be dangerous. 

4.1 Gas density 

Diffuse clouds containing reasonable quantities of molecular 
gas have hydrogen nuclei number densities iih < 100 cm -3 , 
whilst giant molecular clouds have densities ranging from 
an average value of ~200 cm~ 3 up to 10 6 cm~ 3 in their 
cores (see, e.g., Scoville et al. 1987; Glisten 1989; Dyson & 
Williams 1997, and references therein). Therefore, we con- 
sider the sensitivity of X to variations in density within the 
range 10 2 ^ nn ^ 10 6 cm -3 . For the purposes of compari- 
son, we have adopted the simplifying assumption of constant 
density across the PDR in these models. The variation of X 
with visual extinction for densities from 10 2 to 10 4 cm -3 is 
shown in Fig. 3. 

As the density increases from the lowest values consid- 
ered, the minimum value of X drops and the profile min- 
imum gets sharper and shifts to lower Av- The effective 
shielding of CO from photodissociation increases with gas 
density and allows its abundance to reach saturation at lower 
Av values, resulting in stronger emission and therefore re- 
ducing the value of X, which is inversely proportional to 
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Figure 2. Top: an X factor profile calculated for the standard 
parameter values listed in Table 2. CO(1-0) emissivity (A; dotted 
line) and X (solid line) are plotted as a function of Av into the 
cloud. Bottom: CO(1-0) integrated intensity (dotted line) and H2 
column density (solid line) as a function of Av into the cloud. 
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Figure 3. X versus Av profiles for varying number density of 
hydrogen nuclei within the cloud (njj cm -3 ). 



CO emission strength. At higher densities, the CO (1-0) line 
rapidly becomes optically thick with increasing depth, caus- 
ing emission to occur within a thin region of the cloud and 
narrowing the minimum of the X versus Av profile. At den- 
sities above 10 4 cm~ 3 (not shown) the sharpening and deep- 
ening of the X profile minimum becomes more exaggerated. 



© 2006 RAS, MNRAS 000, 1-9 



6 T. A. Bell et al. 



io 2 



io 21 



X 



IO 1 



- \ 1 \ 1 

- \ \ 

- \ I 
\ \ 
\ \ 


i 1 1 

\ 


, , y , , 
\ 






- u=i 

\ 


\ x=io 

\ 


\ x=io 3 


\ x=io 5 

\ 


- 












; x=o.i 


1 


1 


1 





4 6 
A„ [mag] 



10 



Figure 4. X versus Ay profiles for varying incident FUV radia- 
tion field strength (x Drainc). 
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Figure 5. X versus Ay profiles for varying cosmic ray ionization 
rate (£ 1.3 X 10" 17 s" 1 ). 



4.2 Radiation field strength 

Radiation field strengths incident upon gas clouds can vary 
from a few times the standard interstellar Draine field up to 
X ~ 10 7 Draine in regions of intense star formation (see Hol- 
lenbach & Tielens 1999, and references therein). We adopt 
a range of 0.1 ^ \ ^ 10 5 Draine for consideration in this 
study. Fig. 4 shows X versus Ay profiles for radiation field 
strengths within that range. 

The minimum of the profile increases, broadens and 
moves deeper into the cloud with increasing radiation 
strength. The C + /C/CO transition region occurs at higher 
values of Ay as the radiation penetrates deeper into the 
cloud and the rise in CO abundance is more gradual un- 
der stronger radiation, as shielding becomes less effective. 
This leads to a broader emission feature, causing the wider 
profiles seen in Fig. 4 as the radiation intensity increases. 



4.3 Cosmic ray ionization rate 

The cosmic ray flux is known to vary by over an order of 
magnitude. Magnetic field lines can channel cosmic rays 
away from dense molecular cores; alternatively, the flux of 
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Figure 6. X versus Ay profiles for varying cloud age (t yr). 

particles in starburst regions can be many times higher than 
the canonical rate (Schilke et al. 1993). We therefore con- 
sider a range of 0.1 C 10 3 times our adopted canonical 
rate of 1.3 x 10~ 17 s" 1 . 

Increasing the flux of ionizing particles through the 
cloud produces a complex behaviour in the X versus Ay 
profiles (see Fig. 5). As the ionization rate increases from 
the standard value (see Table 2), CO is initially destroyed 
more effectively through reactions with He + and its abun- 
dance drops, causing an upward shift in the X profile. This 
is the case for the £ = 10 and ( = 10 2 models. However, the 
corresponding increase in the cosmic ray heating rate causes 
a rise in gas temperature, becoming the dominant heating 
mechanism at Ay > 1 mag for £ = 10 2 and at all depths 
for (" = 10 3 . This promotes CO formation and emission, and 
counters, in part, the increased destruction rate. A hump 
emerges in the CO abundance and emission profiles (not 
shown) across the region < Ay < 1.5 mag for ionization 
rates above 10 2 , produced by the strong cosmic ray heat- 
ing and effective [Ol] cooling that peaks at Ay ~ 1.5 mag. 
This results in X profile minima that are lower than in the 
standard model and nearer the cloud surface. Thus, a slight 
rise in ionization rate shifts the X profile minimum upwards, 
whilst a further increase in ionization rate leads to the X 
minimum dropping and moving towards the cloud surface. 
In models with ( below the standard value (not shown), 
there is no appreciable change in the X profile shape. 



4.4 Time-dependence 

The evolution of a cloud consisting initially of atomic species 
produces profiles that change from low, very broad minima 
to the usual, sharper minimum as the chemistry approaches 
equilibrium (see Fig. 6) . The formation of molecular hydro- 
gen represents the significant time-scale in the chemical evo- 
lution, since the formation of CO proceeds fairly rapidly 
once H2 is dominant. At early times, hydrogen is mainly 
in atomic form and the CO abundance is low, reaching a 
value of ~10~ 8 for depths beyond the point at which pho- 
todissociation becomes negligible. The low CO abundance 
and higher gas temperature mean that the CO (1-0) line 
does not become saturated and emission extends deep into 
the cloud, producing profiles with very broad minima. Once 
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Figure 7. X versus Ay profiles for varying metallicity (Z/Zq). 
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Figure 8. X versus Ay profiles for varying turbulent velocity 
within the cloud (i>t ur b kins -1 ). 



CO reaches abundances >10~ 5 (which occurs at an age of 
~1 Myr for the standard parameter values) emission from 
the J = 1 — > line becomes optically thick, producing an 
emission peak at Ay ~ 1.2 mag and the characteristic sharp 
minimum in the X profile. For the standard parameters con- 
sidered in this model (see Table 2) there is no noticeable 
change in the X profile for cloud ages greater than 10 Myr, 
even though the chemistry does not reach steady-state until 
100 Myr. This is because emission from the optically thick 
CO(1-0) line is no longer sensitive to the change in CO 
abundance at depths beyond the emission peak. 



4.6 Turbulent velocity 

The degree of turbulence within a cloud can be parametrized 
by its contribution to the observed random velocities. Tur- 
bulent velocities can be close to zero in quiescent clouds 
and 5 kms -1 or more in particularly active regions (e.g. 
Williams, Blitz & Stark 1995; Brunt & Heyer 2002). We 
consider a turbulent velocity range 0.5 < Wturb =S 3 kms -1 
in this work. 

An increase in the turbulent velocity of the cloud serves 
to decrease the value of X at each depth point whilst retain- 
ing the characteristic profile shape (see Fig. 8) . The stronger 
turbulent broadening of emission lines leads to lower opac- 
ities, causing a rise in CO emission, and a corresponding 
drop in the value of X. 



4.5 Metallicity 

Metallicities a factor of 10 below Solar have been observed in 
Local Group galaxies, including the Small Magellanic Cloud 
(Lequeux et al. 1979), and more distant galaxies are believed 
to have metallicities as low as 10 -2 . We adopt a metallic- 
ity range of 10 -2 < Z/Zq ^ 1 for this study. The depen- 
dence of X upon metallicity has been the subject of previous 
observational studies (e.g. Wilson 1995; Arimoto, Sofue & 
Tsujimoto 1996) that find a similar trend to that shown in 
Fig. 7, namely, that the value of X is generally increased 
when metallicity is less than Solar. This is simply due to 
the reduced abundance of CO that can form, resulting in 
a matching drop in CO(1-0) emission. There are, however, 
changes in the profile shape as the metallicity decreases. 
The reduction in metallicity also inhibits H2 formation due 
to depletion of the grain surface area available for cataly- 
sis and this has implications for the time-scales required to 
reach chemical equilibrium. For the Z — 10~ 2 model, the 
abundance of H2 does not reach steady-state until 1 Gyr, 
a time-scale that is unlikely to be reached under real con- 
ditions, where turbulence and dynamical evolution prevent 
clouds remaining quiescent for such long periods. Low metal- 
licity environments can therefore be far from chemical equi- 
librium, in which case the appropriate value of X can be 
significantly different from that derived under the assump- 
tion of steady-state. 



5 CONCLUSIONS 

This work presents a thorough investigation of the proper- 
ties and behaviour of the X factor under varying physical 
parameters. Such an understanding of the variation in X is 
vital when using the CO-to-H2 method to estimate molec- 
ular gas mass in regions that are markedly different to the 
local interstellar medium. We have demonstrated that the 
conversion factor can vary by over an order of magnitude 
under the influence of some physical parameters, in par- 
ticular the density, age and metallicity of the region. The 
benchmarking of the PDR code provides support for these 
predictions, whilst highlighting the inherent uncertainties in 
theoretical models of this kind. To summarize the trends 
discussed in this work: 

(i) Increasing density causes the minimum value of X to 
drop and shift to lower Ay. 

(ii) Stronger radiation fields slightly raise the X profile 
minimum and push it deeper into the cloud. 

(iii) A slight increase in the cosmic ray ionization rate 
will raise the X minimum, but ionization rates over 100 
times the standard value result in the minimum dropping 
and appearing nearer to the cloud surface. 

(iv) Turbulent velocity variation does not significantly al- 
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ter the shape or position of the X profile minimum; increas- 
ing turbulent velocity produces a slight drop in the X value. 

(v) Young, predominantly atomic regions produce broad, 
low X profile minima; the characteristic X profile shape only 
appears once a cloud has undergone considerable chemical 
evolution, corresponding to cloud ages above f Myr. 

(vi) Decreasing metallicity has a strong effect on the X 
profile, increasing the minimum by over 2 orders of magni- 
tude above the canonical value in the Z — 10~ 2 case. 

(vii) The issue of time-dependence in low metallicity en- 
vironments is particularly important, since the reduced H2 
formation rate means that the chemistry does not reach 
steady-state for cloud ages below fOO Myr. 

Conditions that produce extended CO emission deep 
into the cloud (such as low density, strong radiation fields 
or low metallicity) serve to broaden the X profile minimum, 
and the CO emission reliably traces the H2 column density 
to higher Av in such cases. Since interstellar clouds are likely 
to be clumpy on scales significantly less than the f mag of 
extinction considered in these models, the CO emission in 
each clump is likely to be optically thin for most of the 
parameters considered here. Conditions producing extended 
emission allow for the presence of larger clump sizes before 
the reliability of the X factor as a tracer of total H2 column 
density breaks down. 

In future work, we will derive a set of X factor val- 
ues appropriate for specific extragalactic environments and 
compare them to observational studies of the X factor in 
similar regions. This will allow the accuracy of theoretical 
modelling of X to be evaluated and we will present a table 
of conversion factors that may be used by observers to ob- 
tain more reliable estimates for molecular gas masses, along 
with limitations for their use. 
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